Metamagnetism in one dimensional systems with edge sharing CuO polihedra 
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We study a Heisenberg chain with nearest-neighbor (NN) J\ and next-NN J2 exchange interactions 
with anisotropics Ai and A2 respectively. We investigate by analytical and numerical methods the 
region of parameters for which there is a jump in the magnetization M as a function of magnetic 
field B. Some materials with edge sharing CuO polihedra are candidates to show an abrupt change 
in M(B). 



The magnetization as a function of applied magnetic 
field in several materials Q-j^] shows a discontinuity or 
very rapid increase at a certain field B c . Gerhardt et al. 
have shown that for certain parameters, a magnetization 
jump is also present in the spin-1/2 XXZ chain with NN 
and next-NN exchange coupling (keeping Ai = A2 = A) 



H = Y^iMSfSf +1 + SfSf +1 + A l S?S? +1 ) 

i 

i 

-BJ2 S ^ (!) 

i 

For a metamagnetic transition to occur at very low 
temperatures, the zero-field ground-state energy per site 
E as a function of the magnetization M — J2i Sf/L (L 
is the number of sites), should satisfy two conditions: I) 
d 2 E/dM 2 < in a finite interval of values of M. Then 
one can draw a straight line E'(M) which is tangent to 
E(M) in at least two points (the Maxwell construction, 
see Fig. 1 (a)) in such a way that E(M) > E'(M) for 
Mi < M < M 2 . II) E(M 2 ) > E{Mx). If these two 
conditions are satisfied, M jumps from Mi to M 2 at the 
critical field B c = [E{M 2 ) - E{M X )]/(M 2 - Mi). 

From the general behavior of E(M), Gerhardt et al. 
have found that when metamagnetism exists, M 2 = 
1/2 and the condition II ceases to be satisfied when 
Ml = 0. More precisely, from their finite-size results 
for E(M, a, A), with a — J 2 /Ji, they obtained a criti- 
cal value of A (A/(a)) from the equation £7(0, a, Af) = 
E(l/2, a, A/). For A < A/ the system is ferromagnetic 
at B = 0. Another critical value A Q (a) was obtained 



from the condition d E/dM 



M=l/2 



0. For A > A a 



the curvature d E/dM is positive for all M. The dis- 



cretized d 2 E/dM z 



M=l/2 



= has some finite-size effects 



H . From the numerical solution of the problem of two 
spin excitations on the ferromagnetic state for L — > 00, 
more accurate values of A a (cv) were obtained recently 
for a < 1/2 §]. In the region of the (a, A) plane where 
Af(a) < A < A a (cv) a metamagnetic transition occurs 
in the model ||||. 

We have studied the two-magnon problem for generic 



values of Ai and A2, and found analytical results for 
the condition d 2 E/8M 2 \ M=1/2 = if a < 0.75. When 
Ai = A2 = A, in the region a < 1/4, the function Ay (a) 
can be accurately approximated by: 



A f = -l + 2j2 al +6a 5 + 0(a 6 ), if a < 0.2 
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A f = -(-5 + Vl7) - 0.462V1 - 4a, 0.2 < a < -. 
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For 1/4 < a < 1/2, although the algebra is more in- 
volved, the exact result is simpler: 
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Finally in the region 1/2 < a < 0.75, Af (a) is very 
flat. Near a = 1/2 it can be approximated as A/ = 
— \ + 0.309(:r — i) 2 . These results show that metamag- 
netism is not possible if A > (-5 + \/Tf)/4 = -0.219. 
Unfortunately, such a large anisotropy of J 2 seems unre- 
alistic. Instead, Ai = — 1 corresponds to isotropic ferro- 
magnetic Ji, since a rotation of every second spin in tt 
around the z axis changes the sign of the x and y com- 
ponents of J\. 

The main purpose of this work is to extend the pre- 
vious results to negative Ai and positive A2. Since it 
is expected that the parameters for several copper oxides 
containing edge sharing Cu-0 chains lie near the isotropic 
limit Ai = — 1, A2 = 1, H we consider this limit in what 
follows. From numerical diagonalization of 20 sites, we 
obtain that spontaneous ferromagnetism does not take 
place for a > 1/4. If in addition a < 0.7, there is a 
bound state in the two-magnon problem at wave vector 
q 2 = 2qi, where qi = ± arccos[— l/(4a)] are the wave 
vectors of the one-magnon states of lowest energy. For 
a > 0.7, there might be a two-magnon bound state with 
q 2 7^ 2qi, but we have not studied this alternative be- 
cause it seems not possible to solve the problem analyt- 
ically for large a. Thus, we expect a jump in M(B) for 
1/4 < a < a c with a c > 0.7. 

In Fig. 1(a) we show E(M) for a chain of L = 20 
sites with periodic boundary conditions for a = 0.425, 
chosen in such a way that q\ = ±77r/10 are allowed wave 
vectors of the finite chain. For other values of a, one 
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might obtain a numerical negative d 2 E / dM 2 \M=i/2 be- 
cause of frustration effects which increase E(M — 1/L). 
In spite of this precaution, the results show a significant 
even-odd effect: the energies for odd (even) total spin 
S = \M\L seem to be shifted to higher (lower) energies. 
If this effect persists in the thermodynamic limit (keep- 
ing L even) states with odd S become irrelevant (because 
they do not minimize E — MB for any B) and a bound 
state in the two-magnon problem does not necessarily 
imply d 2 E/dM 2 \ M=1/2 < 0. From E(S/L) for the three 
highest even S with L — 28, minimized with respect to 
the optimum twisted boundary conditions to allow for 
incommensurate wave vectors |7j , we obtain a very small 
curvature which is negative for a < a c = 0.359 but pos- 
itive for a > a c - If a c remains finite in the thermody- 
namic limit, M(B) would increase abruptly for a > a c , 
but without showing a true jump. While this difference is 
hard to distinguish experimentally, it would be of interest 
to calculate d 2 E/dM 2 \ M=1 / 2 using larger clusters. 

To obtain a continuous curve E(M) from which M(B) 
can be derived, we have fitted the eleven numerical values 
represented in Fig. 1(a) by a polynomial of even powers 
of M up to M 10 . This function satisfies the physical con- 
dition E(M) — E(—M) and has six fitting parameters 
(nearly half of the number of points to be fitted, to aver- 
age the even-odd effect). The resulting B = dE/dM is 
represented in Fig. 1(b). At a critical field B c — 0.192Ji, 
the magnetization jumps from Mi — 0.347 to M 2 = 1/2. 
While the numerical values of B c and particularly Mi 
depend on the particular fitting procedure used and the 
size of the system, the general shape of M(B) is robust: 
around B/Ji = 0.19±0.01 there is a sudden increase of 
M from - 0.25 to 1/2. While the shape of M(B) does 
not depend very much on a > 1/2, B c increases strongly 
with a. 

To conclude, while the existence of a real jump in 
M(B) requires a study of larger clusters, we have shown 
that the magnetization of the model for parameters ap- 
propriate to edge-sharing Cu-0 chains has a sudden in- 
crease at a magnetic field B c . Using parameters calcu- 
lated for La6CagCu2404i, and assuming a gyromag- 
netic factor g = 2 we obtain B c ~ 14 Tesla. 
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FIG. 1. Energy per site as a function of total spin per site 
for a chain of 20 sites (solid circles). The full line is a fit (see 
text). Dashed line and diamonds correspond to the Maxwell 
construction, (b) Magnetic field as a function of the magne- 
tization. Parameters are J\ = 1, J 2 = 0.425, Ai = — f and 
A 2 = f . 
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